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ABSTRACT 

We study abelian gauge theories with anisotropic couphngs inA + D dimensions. 
A layered phase is present, in the absence as well as in the presence of fermions. 
A line of second order transitions separates the layered from the Coulomb phase, if 
D <3. 



1 Introduction 



The phase structure of gauge theories with isotropic couphngs is known since the 
advent of lattice gauge theories. For U{1) theories, there are two phases: a strong 
couphng and a weak couphng (Coulomb) phase, above and in four spacetime dimen- 
sions, while below this limit there is only one, strong coupling (confining) phase. 

It was realized some time ago []1| that lattice gauge theories with anisotropic 
couplings may have an extra, "layered", phase. The prototype of such a system is 
the 4+ 1-dimensional U{1) theory, in which the plaquettes in four dimensions have 
the same coupling constant (3 = l/g"^, whereas all plaquettes with links in the fifth 
dimension have coupling (3' = l/g'"^- 

In the layered phase all four-dimensional hyperplanes, along the fifth direction, 
decouple from one other and in each hyperplane we have a four-dimensional Coulomb 
phase. This takes place for f3' small enough (typically < 1/d, where d is the dimen- 
sionality of the layer, in our case d = 4) and for f3 large enough (typically f] > 0(1)). 

The intuition of Fu and Nielsen for the occurrence of the layered phase is 
worthwhile bringing up at this point. 

Let us suppose that we are in the 5d Coulomb phase (so /? ~ large). Between 
two test charges we have a Coulomb force in this phase. Imagine making (3' smaller 
but keeping f3 fixed. Then the force between the charges, in the fifth direction, will 
increase and will become, for small enough jS', confining and create a string- whereas 
the force will still be Coulomb-like in the four other directions. Thus we expect the 
Wilson loops to behave like 

• Wf^iy{N^ M) ^ exp — A^Mcr (strong coupling); 

• Wuf,{N, M) exp -(A^ + M)t (Coulomb phase); 1 < /x, z/ < 5 and 

• W^,{N, M) = exp -{N + M)t', 1 < /i, z/ < 4; 
Wf,5{N, M) = exp -NMa' (layered phase). 

There isn't any layered phase with the roles of /3 and f3' reversed. At small enough 
f3 one has either strong coupling or weak coupling and Monte Carlo simulations [Q] 
seem to corroborate this picture, though only single plaquette averages were calcu- 
lated. This asymmetry is due to the asymmetric role of (3 and jS': whereas (3 is the 
coupling restricted to the 4c? space, /5' not only acts in the space orthogonal to the 
Ad space, but also couples the links in the fifth direction to those in the others. 

The layered phase is due to the simultaneous presence of Coulomb forces in the 
physical direction and confining forces in the other directions. So we need at least 
4 dimensional layers for an Abelian theory(i.e. 4 + D spacetime dimensions) and at 
least 5 dimensions for a non- Abelian theory {5 + D spacetime dimensions). 

In this paper we will concentrate on two issues: in section 2 on the order of 
the various phase transitions and in section 3 on the presence of fermions, both by 
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mean field techniques and Monte Carlo methods (section 4). In section 5 we present 
results from mean field and Monte Carlo methods and conclusions and outlook are 
set forth in section 6. 

2 Mean Field Theory for 4+D Abelian gauge 
theory and order of the transition 

Mean field techniques for gauge theories @] typically give first order transitions-in 
contrast to spin systems. The reason may be found in the structure of the fixed- 
point equation for the order parameter, as will be made explicit below and may 
be summarized by the remark that a plaquette has four links (gauge system) while 
there are two spin variables per link (spin system). 

The standard way to do mean field theory in a 4 + D anisotropic system with 
gauge action 

Sg = P Yl {l-m^u{x)) 

l</i<i/<4,x 

a;,l<M<4, 
5<!/<4+D 

+ f3' J2 il-m,,ix)) (1) 

5<fi<iy<'i+D,x 

goes as follows (we use standard notation-the link variables U^j,{x) = e^'^>^^^^ in the 
X ^ X + fi direction, multiplied around a plaquette give the product U^y{x)). 
The partition function Z is 

Z= jX{U{l)e^V-SG (2) 

and may be trivially rewritten, by inserting the identity operator 

1 = / n Dvr{l)Dvi{l) 5{vr{l) - mil)) 6{v,{l) - QU{1)) (3) 
I 

and using the integral representation of the delta function, (5(f—[/) = / daexp —{ia{v— 
U)). The result is 

Z = fY[DU{l)Dar{l)Dai{l)Dvr{l)Dvi{l) 

I 

X exp{iaril)m{l) + ai{l)'^U{l)) 

X exp -iS{v) + iar{l)vr{l) + iai{l)vi{l)) (4) 
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S{v) is obtained by substituting Vr{l) for 3?f/(/) and Vi{l) for 53f/(/) in eq. (0). 

The U{1) integrations now decouple, at the expense of introducing the new hnk 
variables a{l), and give a contribution oc I]j log (-^o {\J + «K0)) the expo- 
nential in eq. (H), leading to an effective action 



5eff = ^(t;) + E log (^0 (^V"'(0 + «K0 j j (5) 

Traditionally one introduces axial gauge in the 4 direction0, U4^{x) = 1 for all 
points. 

The equations of motion for the fields v{l) and a{l) may now be explicitly written 
down. The correct solution is the configuration that minimizes the free energy, 
— logZ. Fu and Nielsen guessed that ai{l) = Vi(l) = for all links / and, 
furthermore, Vf^{x) = v for fi = 1, . . . , 3, while v^{x) = v' for /i = 5, . . . , 4 + D, 
where v, v' are independent of x. With this simplifying Ansatz the equations of 
motions become 

4(3v^ + 2Df3'v% + 2(3v = a 

(6) 

V = Ii{a)/Io{a) 
6(3'v'^v' + 2(3'{D - l)v'^ + 2l3'v' = a' 

(7) 

v' = h{a')/Io{a') 

and the free energy/site (for the D = 1 case, on which we focus in the remainder of 
this section) may be written as 

/ = 3(-i3v^(v^ + i) + av- log Jo(a)) + (3'{-3v%'^ - v'^) + a'v' - log/o(a') (8) 

and delivers plaquette expectation values through 

d ^ 

'oof = X! (^A''^) 

I (9) 

In ref. these equations were analyzed and three regimes identified: 

• f = 0, f ' = (strong coupling, I); 

• f 7^ 0, t>' = (layered phase, II); 



^This is a gauge choice that aUows one to compute the free energy — log Z easily. Formally, axial 
gauge would be Ui+D{x) = 1 for all points x-but the mean field equations become less transparent. 
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V 0, f ' 7^ (weak coupling,!!!). 



The transition on the axis /3 = is an artifact of mean field theory-along this line 
the theory splits into a sum of one-dimensional spin theories. Mean field theory is 
known to break down in low dimensions and corrections shift the "transition point" 
significantly to the right Jll. 

Let us now argue the order of the transition(s) from mean field theory. 

The transitions separating ! and !!, ! and !!! are first order. This is the expected 
result from mean field theory, since the term in eq. (|^) leads to a jump from v = 
(!) to f 7^ (!!,!!!). For the strong to layered transition, our result is in agreement 
with the latest Monte Carlo data, obtained for /?' = 0, i.e. four- dimensional compact 
f/(l) theory 0. Along the (!!,!!!) line, however, v stays non-zero; the transition is 
marked by v' becoming non-zero. The mean field equation (|^) is now relevant -and 
for D = 1 v' enters linearly, due to the fact that f/5 comes in only quadratically in the 
4+1 dimensional anisotropic model! The transition is second order-a.nd Monte Carlo 
results 1^, ^ do indeed confirm this. For D > 1 this argument is no longer valid, 
since v' enters cubically-however, the coefficient of this term is not large enough to 
drive the transition first order-at least as long as -D < 3. Monte Carlo results are 
not yet available for the higher dimensional case {D > 1), so we cannot check the 
validity of this mean field prediction. 

3 Coupling fermions 

The presence of fermions renders the order parameter inefficient -it will always indi- 
cate perimeter behavior. What one should look at instead is the fermion spectrum, 
which is different in the three phases: 

• bound states with strong coupling (!); 

• fermions move freely in the four-dimensional subspace as in a four- dimensional 
Coulomb phase, but feel a confining force in the extra directions (!!); 

• free fermions, positronium states in A + D dimensions. 

We have checked, by mean field methods, that the phase diagram described above 
is not qualitatively altered by the presence of the fermions. In this section we will 
explain these methods. 

Our starting point is the fermionic action, with a Wilson term 



Sf 



- ^ ^(x)7^[f/^(x)^(x + fi)- f/t(x - fi)^{x - fi)] 



+ 
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X 

= (10) 

x,x' 

The total action is the sum of the above and the gauge action, equation (|l|) 

S = Sg + Sf (11) 

The partition function now reads 

Z = [llD'^{x)D^{x)l[DU{l)exp-S (12) 

X I 

The fermions may then be integrated out, leading to the effective action 

S,fi = S-TT\ogM{U) (13) 

The mean field approximation goes through exactly the same steps as in the pure 
gauge case, leading to 

/ioo r 
l[Da{l) Dv{l)exp-{Ses + (yv -\oglo{a)) (14) 
-ioo ^ J 

The difference with the pure gauge case is that the fermionic determinant leads to 
a non-local action for the f-variables. The equations of motions now become 



ipv^ + 2Dp'v% + 2(3v + ]{v, v') = a 
V = /i(a)//o(a) 

and 

Gp'v^v' + 2(3' {D - l)t7'3 + 2(3'v' + ]'{v, v') = a' 

v' = h{a')/Io{a') 
The currents j{v,v') and j'{v,v') are short-hand for 



(15) 



(16) 



jiv,v') = ^_^^^ TrlogA<|„^(^)=„,„,(:,.)=^/, /i = 1,2,3; z/ = 5,...,4 + L' (17) 

and 

j\v,v') = j—--Tr\ogMl^^^)=^,^,^(^)=^', yu = 1,2,3; z/ = 5,...,4 + D (18) 



So j{v,v'), j'{v,v') are the expectation values of the components of the fermionic 
current in the background field {v,v'). Straightforward manipulations lead to the 
expressions 



Tr (^G 



V, V 



+ Tr 



f , V 



/i= 1,2,3, 5,...,4 + D 



(19) 



where G{v,v') is the fermion propagator in the background {v,v') and the trace is 
only over Dirac indices. 

Since the background is constant in space, the currents can be computed 
implicitly as functions of v, v' through a 4 + D dimensional momentum integral. 

Let us now introduce the variables 



and 



/ 3 4+D 

PF = M — rl^(l — vcospx) + 1 — cosp4 + X! (-'^ ~ v' cospx] 

\A=1 A=5 



3 4+D 

P = '^v'^ sin^ Px + sin^ ^4 + ^ f sin^ px + W"^ 

A=l A=5 



in terms of which the currents may be written as 



J, = 4 



and 



3u 



. (27r)4+^ 



fsin p^ + r cosp^W 



P 



/i = l,2,3 



(27r) 



4+D 



v' sin py + r cos pi,W 



P 



z/ = 5, . . . , 4 + D 



(20) 



(21) 



(22) 



(23) 



Obviously, if w ' = then j' = (c.f. equation ([T6|) -or, more explicitly, by combining 
(pop and (p3|)) and the equations of motion are identical with or without fermions. 
Assuming a second order transition between the layered and the Coulomb phases 
we find the shift of the transition point induced by the presence of the fermions (for 

f = 1) as 

Pc=\-^ (24) 



where 



TTJn (27r)^ 



4 8 

p2 



(25) 



where S = Z]^=i sin^p^, W and P as above, but for t>' = 0,f = 1. Note that the 
shift is negative and very small (except for special values of M where P has poles). 
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Figure 1: The current j, as a function of for (3 = 1.2. 



i.e. the layered phase, in the presence of fermions, loses ground against the Coulomb 
phase. Physically the sign may be understood from the fact that the presence of 
fermions weakens the confining forces, responsible for the layered phasdl. Thus there 
is a critical number of fermions, beyond which there isn't any layered phase. The 
fermionic currents j are independent of (3' in the layered phase (since v is). The 
current j' is identically zero in that phase. For fixed (3 it grows in the Coulomb 
phase, cf. figs. |l], 0. 



4 Monte Carlo Methods 

The size of the five dimensional system was taken to be 4^ and 8^. For the pure 
gauge simulations we used a 5-hit Metropolis algorithm with dynamically adjustable 
stepsizes e,p and in order to maintain an acceptence rate of roughly 50% in each 
direction. 

As we are interested in the phase structure, we performed three thermal runs on 
each lattice size. Such a run consisted of 5000 updates at the starting values of the 
couplings P and for thermalization, followed by 1000 updates and measurements, 
after which either (3 or (3' was increased by 0.01, and another 1000 updates and 
measurements were made without rethermalization. After having reached the final 
values of (3 and (3', the process is reversed intil the starting values have been reached 

"in our original paper |^ there is an unfortunate error in the drawing of the phase diagram. 
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Figure 2: The current js, as a function of /?', for /? = 1.2. 

again. There were no thermahzation sweeps done at the start of the second leg of 
the run. 

The number of 1000 measurements was chosen with the following reasoning in 
mind. In case of a first order phase transition, there are, loosely speaking, two 
autocorrelation times at play: the first being the tunnelling time r^, and the second 
being the relaxation times within the states i. The number of measurements, 
N^eas, should then satisfy 

< N «: 

' l \ -i ' meas ' 

For second and higher order phase transitions, t'^ is absent, and we are left with the 
condition 

< N . 

' 1 -L 1 mcas • 

Although changing the couplings will bring the system out of equilibrium, we 
might hope that this disturbance is relatively small, and the system will relax to 
equilibrium within sweeps. 

As observables to measure we choose the average plaquettes U^p and C/5, 

t^sp = E ^U,^i^)) (26) 

IJt<V<5 
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C/5 = 



{4^ E ^u^(-)) 



(27) 



At<5 



and the Polyakov line correlators P^p and P5 




(28) 



//<i/<5 




(29) 



where Pij,{x) is the Polyakov line starting at point x, running in the ft direction and 
N is the lattice size in the corresponding direction. As the phase transitions in this 
model are conjectured to be confining-deconfining transitions in one or more dimen- 
sions, these observables are indeed order parameters for these phase transitions. For 
a confining phase, P will be zero as the system will be disordered over long distances. 
In a Coulombic phase, however, P will be non-zero, as now such long correlations 
are allowed. 

When including fermions, we used a Hybrid Monte Carlo algorithm with a first 
order discretization of Hamilton's equations. A standard Conjugate Gradient algo- 
rithm was used for matrix inversion. One sweep consisted of momentum refreshment, 
followed by 10 integrations with stepsize r = 0.05. An accept/reject step was done 
in order to maintain balance. The runs were performed as in the pure gauge case, 
but now with 1500 thermalization sweeps and 250 sweeps at each (/?, /?') pair. 

The actual simulations were done on the Wuppertal CM-5 in its various incar- 
nations. We are, therefore, unable to make comments about the performance of the 
two programs on this machine. 

5 Results 

This section will be subdivided into two parts: one for the pure gauge case and one 
including the fermions. For each case we shall confront results (namely plaquette 
values) from mean field theory with Monte Carlo simulations-so we work in five 
spacetime dimensions in what follows and will comment on what happens for D — 
2, 3 at the end. 

5.1 Pure gauge results 

In this paragraph we report on mean field calculations and Monte Carlo simulations 
for the pure gauge case. The quantities of interest are the average plaquettes. 
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Figure 3: Average space-space plaquette vs. /3 at = 0.2. Solid line is mean field 
theory, points are Monte Carlo results from 1000 measurements on an 8^ lattice. 

Using the equations of motion, (|]) and (0), one obtains the following expressions, 
within mean field theory 

ld_ v' + v' 
, _ I d ^_ v'\3v^ + 1) ^^"^ 

These reflect the fact that the equations of motion cancel all implicit dependence 
of the free energy on f3 and So the rhs of eq. (|30D represents only the explicit 
dependence and is, therefore, equally valid when fermions are included. In fig. |^ we 
display typical plaquette values vs. P at fixed relevant for the strong ^ layered 
transition. Mean field theory predicts a transition at /? = 0.8, while the simulations 
indicate that the transition is at /3 = 1.02. It should be noted that corrections 
to mean field theory also indicate an upward shift in the transition point In 
the layered phase, the mean field prediction is that the space-space plaquette, for 
fixed /?, doesn't depend on /?'. The reason is that in the layered phase f ' = and 
/?' enters the mean field equations only in conjunction with it; furthermore, from 
eq. (|30|), the plaquette in the fifth direction is zero. There is, however, a strong 
coupling correction 6U = P'/2 + 0(/?'^) that must be added-this is a well-known 
result that may also serve as a check on the simulations. 

Finally, we present evidence that the layered<-i> Coulomb is second order, namely 
hysteresis loops, in fig. 11. As a comparison, we display hysteresis loops for the strong 
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0.28 0.282 0.284 0.286 0.288 0.29 0.292 0.294 0.296 0.298 0.3 



Figure 4: Hysteresis loop for the layered ^ Coulomb transition at /3 = 1.2. Mean 
field theory. 

^ Coulomb transition in fig. |^. 

For the strong<->layered transition we also display spatial Polyakov line data in 
Figs. ^. These data are the first Monte Carlo data using the order parameter, not 
just plaquette expectation values, that show the occurrence of the layered phase. 
Our results may be summarized by the phase diagram for the pure gauge theory, 
displayed in fig. |^. The main novelty is the evidence that the layered ^ Coulomb 
transition is second order, in agreement with the intuitive argument of the previous 
section. For higher dimensions of spacetime, relevant for other theories, the situation 
is the following: mean field theory predicts that the layered Coulomb phase 
transition stays second order for A+D < 7, while it becomes first order for 4+D > 7; 
this is based on a calculation of the hysteresis loops. Furthermore, as may be seen 
from fig. ^, the strong coupling phase, for 4 + D = 6 is significantly reduced, in favor 
of the Coulomb phase. Unfortunately, Monte Carlo data are not available to check 
whether this is an artifact of the mean field approximation. 

Of course there are several caveats that must be kept in mind: (a) corrections to 
mean field theory may be important-the analysis of ref. [|I| indicates that they do not 
lead to qualitative changes in five dimensions, but are significant in three dimensions; 
regarding these corrections, it should be noted that they tend to soften the first 
order transitions usually encountered (b) regarding the numerical simulations, 
the analysis of ref. 0of the order of the transition is certainly relevant here-the 
results presented are preliminary, certainly indicative of an issue that requires more 
extensive work |^ and, of course, a very long correlation length cannot be ruled out. 
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0.6 



Figure 5: Hysteresis loop for the strong ^ Coulomb phase transition at /3' = 1. 

5.2 Fermions 

We use the Wilson action, as discussed in the previous section, with one fermion 
flavor, and set r = 1 and M = 1. The entire analysis goes through intact-the only 
calculational complication arising from the fact that every iteration of the mean field 
equations towards the fixed point requires the numerical evaluation of a four (in the 
layered phase) or five dimensional integral for which we used the Vegas routine [§. 
Computing, once more, the hysteresis loops through the transitions, we find that the 
presence of the fermions doesn't change the order of the three transitions. Consistent 
with this result we find a shift of the layered ^ Coulomb transition in favor of the 
Coulomb phase. 

To check these predictions we have performed Monte Carlo simulations, measur- 
ing the same observables as in the pure gauge case, using the Hybrid Monte Carlo 
method and our results are summarized in the phase diagram in fig. The caveats 
stated above apply here as well. High statistics simulations are needed to clarify 
whether the results obtained for the layered Coulomb transition do indeed hold 
up. The claim is that the problem deserves this effort 0. 

5.3 Continuum Limit 

The presence of the line of second order transitions in the phase diagrams is in- 
triguing. On both sides of this line we have a phase with zero mass particles and a 
Coulomb law; the difference is that the force falls off like in the four physical 
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0.85 0.9 0.95 1 1.05 1.1 1.15 

Figure 6: Spatial Polyakov line correlators at (3' — 0.2 from an 8^ lattice 
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Figure 7: Phase Diagram of the five dimensional, U{1) pure gauge theory. 
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Figure 8: Phase Diagram of the six dimensional, U{1) pure gauge theory. 
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Figure 9: Phase diagram including fermions of mass M = 1 (points). Mean field 
theory predicts the strong to layered transition at /3 = 0.8; Monte Carlo estimates 
it at 1.02. The pure gauge phase diagram is included for comparison (solid line). 

directions in the layered phase and like l/r^"*"^ in all A + D directions in the weak 
coupling phase. 

In the layered phase the Wilson loops in the extra directions obey an area law. 
At (3' = the string tension becomes infinite and no correlation in these directions 
survives-we have a stack of decoupled Ad layers, each of unit thickness (for D = 1 
or "volumes" for D = 2,3). At /3' > , but /3' < (so v' = 0), one has finite, 

but exponentially damped, correlations between the layers. 

Presumably a continuum limit of this system exists at the tricritical point in 
figs. 0,PJ^. At P' = we cannot take any continuum limit 

6 Conclusions and outlook 

We have presented evidence, analytical (within the mean field approximation) and 
numerical (through Monte Carlo simulations) that gauge theories, with anisotropic 
couplings, naturally support layered phases that admit a non trivial continuum limit. 

Our results show unambiguously a line of second order phase transitions, be- 
tween layered and Coulomb phase in fig. |^. Furthermore, this statement remains 
true in the presence of fermions, fig. ^ The implications of these findings are quite 
broad in scope. In the same context one may find a natural setting for investiga- 
tions of membranes and models with chiral fermions. For this last case the present 
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calculations are a necessary first step for the location of the chiral layer [Q, 

In this paper we only looked at abelian theories, in 4 + D dimensions. If we heat 
up the system, one can include non-abelian systems, as soon as they drop into the 
deconfined phase. We would expect our results for the Polyakov loops to remain 
qualitatively the same-with the caveat that only the Polyakov loop in the Euclidean 
time direction deconfines. Systems of this kind have been studied by other means 



in ref . |[T0|| 
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